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We study the collective modes in relativistic electromagnetic or quark-gluon plasmas with an 
asymmetry between left- and right-handed chiral fermions, based on the recently formulated kinetic 
theory with Berry curvature corrections. We find that there exists an unstable mode, signaling the 
presence of a plasma instability. We argue the fate of this "chiral plasma instability" including the 
effect of collisions, and briefly discuss its relevance in heavy ion collisions and compact stars. 
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Introduction. — Parity-violating effects related to the 
quantum anomalies play important roles in a wide range 
of physics from quantum Hall systems to the cosmol- 
ogy. One example in the transport phenomena is the 
parity-violating current in the presence of a magnetic 
field and an asymmetry between left and right-handed 
fermions, characterized by the chiral chemical potential 
/is = /iij — fiL. This is called the chiral magnetic effect 
(CME) [1-4]. Recently hydrodynamics [5] and kinetic 
theory [6-8] have been appropriately modified to describe 
quantum anomalies and the CME (see also Refs. [9, 10]); 
in the kinetic theory, essential corrections are the Berry 
curvature, the concept diversely applied in condensed 
matter physics [11]. These developments enable us to 
systematically understand the properties and dynami- 
cal evolution of various systems with hitherto neglected 
anomalous effects. 

In this paper, we study the collective modes and their 
consequences in relativistic electromagnetic or quark- 
gluon plasmas (QGP) at finite temperature T and ^5 
based on the new kinetic theory. (For a review on the 
collective modes without Berry curvature corrections, see 
Ref. [12].) We show that the Berry curvature corrections 
dramatically modify the dispersion relation of the collec- 
tive modes, and in particular, lead to unstable modes, 
signaling the presence of a plasma instability. We shall 
call it the chiral plasma instability. Closely related insta- 
bility was studied (but not by the kinetic theory) within 
the electroweak theory at large lepton chemical poten- 
tial [13, 14] and in the context of the early Universe at 
T ^ /15 [15-17] where it is ascribed to a possible origin 
of the primordial magnetic field. 

Our main purpose of this paper is to reveal the poten- 



tial importance of the chiral plasma instabilities in heavy 
ion collisions and compact stars. It has been argued that 
the QGP created in noncentral heavy ion collisions may 
contain locally finite fi^ for quarks and a large external 
magnetic field to yield an observable CME [18]. It was 
also suggested that a degenerate electromagnetic plasma 
with finite for electrons may exist inside neutron stars, 
where the CME is generated in a strong magnetic field 
[28]. We show within the kinetic theory that these QCD 
and QED plasmas at finite fi^ are dynamically unstable 
and reduce ^5 by converting it to (color) electromagnetic 
fields with magnetic helicities. We estimate the typical 
time scales of the chiral plasma instabilities for the QCD 
and QED plasmas as tqcd ^ 1/(q^sM5 Ina^^) for /is ^ T 
and tqed ^ 1/(q^^M5)j respectively. 

Kinetic theory with parity-violating effects. — Let us de- 
scribe a chiral plasma within the kinetic framework. We 
first consider the collisionless kinetic theory and we shall 
argue the effect of collisions later. We also consider 
the regime of a sufficiently weak gauge field A^^, where 
there is no essential difference between Abelian and non- 
Abelian gauge fields. For simplicity of notation, we first 
consider QED, and we will generalize it to QCD later 
(with some modifications mentioned below). We assume 
massless quarks and spatially homogeneous /Z5 below. 

Recall the Maxwell equation 
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where jj'^^ is the induced current and j^^t is the external 
current. For the small gauge field A'^, the induced current 
can be expressed, via the linear response theory, as 



jUK)=Ii'''iK)A,{K), 



(2) 
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in the momentum space, where 11^'' is the (retarded) self- 
energy and K^^ = (w,k) is the four-momentum. From 
Eqs. (1) and (2), we obtain 



Jext ■ 



(3) 



Due to the gauge invariance, which allows us to shift 
+ Kjy, we can choose the temporal gauge = 
0. Then Eq. (3) can be rewritten, using the physical 
electric field, as 



(4) 



where k = |k|. The dispersion relation for the collective 
modes in the system can be found by computing the poles 
of A*-'; a mode with the dispersion relation lo = io{k) 
satisfying Im(u;) > 0, if exists, implies an instability. 

In order to find the explicit form of 11^'^ including the 
anomalous parity- violating effects, we use the kinetic the- 
ory with Berry curvature corrections [6-8] . Although the 
same result for 11^'^ may be obtained in perturbation the- 
ory for the isotropic equilibrium distribution function rip 
[8] (see also Ref. [19]), the kinetic framework would be 
more versatile; it is easy to include the effect of collisions, 
as we shall do later. It is also suitable for other future 
applications, such as the inclusion of the anisotropy of 
rip and numerical simulations of the dynamical evolution 
of plasmas. 

At the leading order in A^, the kinetic theory is given 

by [8] 

{dt+v Vx) np + {eE + vxeB- V^Cp) • VpUp = 0, 

(5) 

where v = p/p, Cp = p{l — eB • ftp) with rip = ±p/2p^ 
being the Berry curvature for right- and left-handed 
fermions, respectively. The energy of the chiral fermion 
is shifted from Cp = p by the amount — peB • flp due 
to the magnetic moment of chiral fermions at finite n [8] . 
This correction makes the term involving VxCp in Eq. (5) 
nonvanishing for an inhomogeneous magnetic field, unlike 
the conventional Vlasov equation. In the presence of the 
Berry curvature flux, the definition of the current is also 
modified to [6-8] 

j = -e^ / [cpVpTip -t- {ftp ■ VpTip) epB -|- Cpftp x Vxrip] . 
Jp 

(6) 

For a moment, we concentrate on right-handed 
fermions with chemical potential /i. (This theory coupled 
to dynamical gauge fields itself is not well-defined due to 



the gauge anomaly, but we shall eventually consider the 
theory with both right- and left-handed fermions, so that 
the gauge anomalies are canceled out.) We solve the lin- 
earized equation of (5) in terms of the deviation Sup from 
the thermal equilibrium state rip = l/[e(^p^^^/^ + 1], 
where Hp = 71° -f Sup. Note that rip does not obey 
rip = n° p and breaks parity. 

Substituting the solution rip to the current j and us- 
ing the linear response theory (2), the self-energy of the 
gauge field can be expressed as Il^''{K) = H'^(K) + 
nr(/^), where [8] 



ri/O 



dv 
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m{K) = ^ie'^'k'^ ( 1 - - ) [1 - ^L{K)] 
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are the parity-even and parity-odd self energies with 
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(10) 



Here denote the spatial indices {K) is vanish- 

ing otherwise] and v = (l,v). Note that n_ is shown 
to be T-independent from the topological nature of the 
Berry curvature (see the similar computation of the CME 
in Ref. [6]). 

Collective modes. — We perform a tensor decomposi- 
tion for the self-energy 11'^. Recall that 11'^ is not sym- 
metric with respect to the indices i and j due to the 
parity- violating term IIj^; we need to use not only the 
longitudinal and transverse projectors, Pf^ = fc'fc-' and 
(5*-' — k^y, but also the antisymmetric tensor, 
l^ijkj^k^ to fully decompose the tensor structure 
of n*-', where fc' = fc'/fc. Then [A~^]*-' can be written as 



pij 



Cl = 

where 
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[l-wL(K)]. 



(11) 

(12a) 
(12b) 

(12c) 



In order to compute A*-' [the inverse of Eq. (11)], it is 
convenient to use the properties of the (mutually com- 
mutable) projectors. 



Pl 



Pt 



PlPt = 0, 



Pi = Pt, PlPa = 0, PtPa = Pa, 



(13a) 
(13b) 
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FIG. 1. Real and imaginary parts of the dispersion relation 
oj — uj{k) for the unstable mode. 



from which we obtain 



1 



-Pi- 



Ca 



C^n- (14) 



Therefore, the dispersion relations of collective modes are 
given by Cl = and Ct ± Ca = 0, which respectively 
reduce to 



nT{K)±JiA{K). 



(15) 
(16) 



While the dispersion relation for the longitudinal com- 
ponent is identical to the case without II^i, that for the 
transverse component is split up into two parts by the 
effect of Ha ■ 

Let us compute the dispersion relations of the collective 
modes based on Eqs. (15) and (16). 

In the long- wavelength limit |a;| ^ fc, II^i is negligible 
compared with IIl,t in the leading order, and 



n 



1 



L.T 



(17) 



The dispersion relations for the longitudinal and trans- 
verse modes are both given by — Wpj, implying that 
nonstatic gauge fields oscillate with the plasma frequency 

Wpl. 

In the quasi-static limit ^ fc, using ujL(K) = 
T^x + x'^ + is"* + • ■ • for positive and negative Im(w) 
with X = <C 1, the transverse dispersion relation 

can be obtained perturbatively in x. Equation (16) with 
the minus sign has the solutions. 



w = ±z7(fc), 7(fc) = 



4a^^,|^_^ (18) 
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Trfc 



for < fc < a/i/TT with a — e^/(47r), while that with 
the plus sign does not; see Fig. 1 for the numerical result 
of Eq. (16) in this regime. The mode uj = ij(k) with 



7 > is unstable and grows exponentially as e"*"^* ~ e^*, 
indicating the plasma instability. The typical scales of 
the unstable mode are 



k ~ afi, 



(19) 



and the time scale of the plasma instability is Tinst — 

1/7 l/(aV). 

Without Ha originating from Berry curvature cor- 
rections, there was no unstable mode for isotropic rip. 
Note also that the criteria for the conventional (Weibel) 
plasma instability in Ref. [20] is not applicable here. 
This is because the assumptions there, = n^p and 
n'J = n^', are violated by the parity- violating effects. 

Now these results can easily be extended to a plasma at 
finite vector and chiral chemical potentials, fxy = fJ-R+fJ-L 
and /i5 ~ /ifl — fiL, by the following replacements: /x^ 
I^'r + A^ U^^v + in ™Z5> and fi ^ fi^ - fiL ^ fi^ in 
n_, and Eq. (18). 

The fate of chiral plasma instabilities. — What is the 
fate of chiral plasma instabilities above? Below we shall 
argue that this instability reduces /xs , so that it is weaken. 
We first recall the relation of triangle anomalies. 



dtQ5^ — I E- B, 



where is the global chiral charge. Equivalently, 



di 



(Q5 + ^h) =0, J AB, 



(20) 



(21) 



where H is the Chern-Simon number (which is also called 
the magnetic helicity in plasma physics). Equation (21) 
is the conservation of the helicity. 

We now explicitly show the reduction of fi^ within the 
kinetic theory in the regime under consideration. We 
concentrate on the unstable mode, Ct — Ca at |a;| ^ 
k, for which the longitudinal part is negligible. Using 
Eq. (14), = A*Jj4t [see Eq. (4)] reduces to 

il3{k) 



(22) 



UJ — ij(k) 

where we defined (3{k) = —2k/{nm?j-,) and J^^^iK) = 
{Pj! - P^)ji^^{K). In the (t,k) space, satisfies 

[dt~i{k)]A\t,]^)=mJUt,^)- (23) 

This can be solved in terms of A^ for jcxt(^, k) = 
Jcxt(k)<5(t): 

A'(i,k)=/3e^*0(t)J^,t(k). (24) 
We then obtain ^ -dtA' and = -kP'^A^ as 

E\t,k) = ~(3[^e'''0{t)+S{t)]J:M, (25) 
B\t,k)=kPe-''e{t)Jl^,{k), (26) 
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where we used P^Jcxt = ^Jc^t that follows from 
Eq. (13b). Hence we find E(i,k) • B*{t,k) < and 
A{t, k) • B*{t, k) > 0. From Eq. (20), this unstable mode 
decreases /j.^ for fixed T and increases the magnitude of 
"H. As a result, it weakens the instability. 

The concrete time evolution of the plasma beyond the 
leading order in A'^ can be described by the full kinetic 
theory [see Eq. (15) in Ref. [8]] together with Eqs. (1), (6) 
and (20), which woiild require; a numerical analysis. Still 
the saturation of the instability itself may be understood 
from the energy and helicity conservations. We assume 
fi5 ^ T. As the energy and magnetic helicity of electro- 
magnetic fields, ~ (kA)"^ and ^ akA"^, come from those 
of chiral fermions, 0(/x|) and ©(/ui), we get the typical 
scales of k and A relevant to the instability: 

kr^anQ, B^iil. (27) 

The electromagnetic fields cannot grow beyond this order 
and will saturate in the end. At the saturation, the gauge 
fields become nonperturbatively large. 

Effect of collisions. For /is ^ T , the mean free time 
between scattering of the plasma constituents in QED, 
Tmft ~ l/(Q!^rina~^) [21], is shorter than the time scale 
of the plasma instability, Tinst ~ ^/{ct^lJ-b)- On the other 
hand, the mean free path, Zmfp ~ Tmft, is larger than the 
typical wavelength of the plasma instability, l//ci„st ~ 
l/{afi^), and one expects that the effect of collisions is 
irrelevant for it. 

To see this explicitly, consider the effect of collisions 
in the relaxation time approximation: we add the term 
— (Jrip/rrei in the right-hand side of Eq. (5), where Trei 
is the the relaxation time which is assumed to be of the 
same order as Tmft. Then lo in the denominator of Eq. (7) 
is replaced by w+i/rrei- Repeating a similar computation 
to above, one finds that the modification to Eq. (18) by 
the effect of collisions is subleading in l/(A;instTrei) <C 
1 and is negligible. (From this argument, the effect of 
collisions is clearly negligible for fi^ ^ T.) As we shall 
see below, however, this is not the case in QCD. 

Quark-gluon plasma with /is. — Suppose /^s (~ T) is 
initially generated in the QGP before the thermalization, 
whose evolution is described by the kinetic theory. The 
above arguments are then extended to color electromag- 
netic fields, within the leading order in A, by replacing 
a ^ as = .gV(47r) and mf, (TV/ + 2Nc)g^T^/6 + 
Nfg^n'^/{2TT^) with Nf and Nc being the number of 



flavors and colors. The qiialitativc difference from the 
QED plasma is that the mean free time (path) is now 
Tmft ~ ^mfp ~ l/(Q;srinQ;J^) due to the fact that col- 
ored gluons, exchanged between quarks, can change the 
color charges of the scatterers [22-25]. As Tmft •C Tjnst 
and Imfp *S ^/hnst in this case, the collisions can affect 
the plasma instability. In the relaxation time approxima- 
tion, L{ljj + i/Trei, k) ~ —iTrei, and one indeed finds that 
Eq. (18) is modified to 

2iv^/_^y 

The unstable modes in < fc < as ^5 /it still exist but 
grow more rapidly with rinst ~ 1/(cKsA''5 l^aj"'^). 

Conclusion and outlook. — We have shown that the rel- 
ativistic electromagnetic plasma and QGP with finite /X5 
are dynamically unstable and jjL^ is damped by the ex- 
ponential growth of (color) electromagnetic fields at the 
early stage [26]. As the gauge field grows, the nonlin- 
earity of the gauge field becomes important: the detailed 
analysis for the QGP would require the kinetic theory 
with non-Abelian Berry curvatures, which is beyond the 
scope of the present paper. 

To understand the roles of the chiral plasma instabili- 
ties at the later stage more quantitatively, detailed stud- 
ies of the dynamical evolution of plasmas based on the ki- 
netic theory need to be worked out. There are also other 
directions in which one can extend or improve our anal- 
ysis: (i) It is straightforward to include the anisotropy of 
rip analogous to Ref. [29] . This allows for the study of the 
coexistence (or competition) with the Weibel instability, 
(ii) The effect of the quark mass can be incorporated 
in the chirality evolution in Eq. (20), which makes /is 
smaller. These analyses would be important to under- 
stand whether the CME remains an observable effect in 
the realistic situations. 
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